r av 

MOf 




NPS-59NN72081A 



NAVAL POSTGRADUATE SCHOOL 



Monterey, California 



ANALYSIS AND SIMULATION OF A 
SUBMARINE STEERING CONTROL SYSTEM 

by 

Ernest Everett Wessman 
and 




Robert H. Nunn 



1 August 1972 



p Distribution unlimited. (See DTIC) 

f 

g 



Only; 

quests 

ost- 



e 023. 



FEDDOCS 
D 208.14/2: 
NPS-59NN72081A 



DOWNGRADED 

APPROVED FOR PUBLIC RELEASE 




NAVAL POSTGRADUATE SCHOOL 
Monterey, California 



Rear Admiral Mason Freeman M. U. Clauser 

Superintendent Provost 



ABSTRACT: 

An existing computer simulation of a hydraulic submarine steering 
control system was analyzed to determine the source of its unrealistic 
behavior, which included verifying mathematical equations and numerical 
parameters, duplicating the simulation on another computer, and investi- 
gating the validity of the assumptions made in its development. Rudder 
turning rate and load pressure drop oscillations which were thought to 
be unrealistic were caused by the neglect of significant rudder moments 
related to the control surface turning rate. Knowledge of actual system 
performance was found to be incomplete and further research was suggested. 
Potentially beneficial modifications in the servovalve and actuator 
simulators were identified. Continuing work in the field of unsteady 
rudder behavior was suggested in order to improve simulator and actual 
system designs. 
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I. THE NATURE OF THE PROBLEM 



A. BACKGROUND INFORMATION 

In 1971, the Naval Ship Research and Development Center (NSRDC) , in 
Annapolis, Maryland, identified a potential stability problem in an 
electrohydraulic servosystem which was being developed. The control 
system was being designed to control the movement of a submarine's rudder 
in accordance with commands from the helmsman. The NSRDC developed an 
analog computer simulation (hereafter referred to as the "baseline simu- 
lation") which incorporated many of the characteristics of an actual 
hydraulic control system, including friction, hydraulic and mechanical 
characteristics, servovalve behavior, and interdependence of elements 
necessary to create a simulation of the actual system. The simulation 
was based on equations developed to model an existing system used in 
the USS PARGO (SSN 650) . Provisions were made so that the simulation 
could be adjusted by changing values of the various components of the 
model, in order to conform to possible changes of the parameters of 
the system under development, such as changes in dimensions, materials, 
and operating conditions. The provision for adjustment of the simulator 
parameters also allowed simulation of other existing hydraulic systems 
of this type. All parameters which were assumed to be of significant 
impact on the performance of the system were incorporated into the 
simulation. 

When the simulator parameters were arranged to model the performance 
of the system under development, the model did not conform to expected 
system behavior; oscillations in pressures and control surface turning 
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rates were observed which had not been observed in similar systems in 
use. Further investigations failed to reveal any obvious errors in the 
simulation, or deficiencies in the parameters which were believed to be 
of significant impact on the performance of the system. 

This report is a result of a project initiated in support of the 
simulation attempt at the NSRDC. The purpose of this project was to 
determine the reasons for the observed behavior of the hydraulic system 
model, and to develop appropriate modifications and refinements of the 
baseline simulation, in order to more closely model actual system 
performance. 

B. SYSTEM UNDER STUDY 

1 . Introduction 

The following paragraphs describe the general characteristics of 
the idealized system that the NSRDC used in the development of its simula- 
tion. The description does not completely describe any real system. 

Rather, it is a model which was designed to include all of the steering 
system characteristics that were expected to be of significant impact on 
the performance of the system. 

2. General System Descriptions 

The steering system (Figure 1) included the servovalve, connecting 
tubing, actuating ram, mechanical linkages, rudder, and minor components. 
The operation of the system is described below. (Numbers in parentheses 
refer to corresponding labeled parts on Figure 1.) 

A desired change of rudder angle is initiated by moving the control 
wheel (1) by hand, which generates an electrical signal having a voltage 
proportional to the degree of movement. This voltage representing the 
desired control surface angle is then compared with a voltage proportional 
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RETURN 

FIGURE I 

SUBMARINE STEERING CONTROL SYSTEM 



to the. actual position of the hydraulic ram (6). An error signal is 
generated which approximates the difference between the desired and 
actual rudder positions (2). This signal actuates an electric torque 
motor (1) in the servovalve, which moves a small flapper (4) located in 
the narrow space between two nozzles, causing an increase in the pilot 
pressure on one nozzle, and a decrease on the other. This difference 
in pressure causes the spool valve (5) to move, thereby permitting 
pressure and subsequent fluid flow to be applied to the hydraulic ram 
(6), resulting in a correcting movement of the control surface (7) through 
the mechanical linkage. The value of the error signal is modified by 
movement of the hydraulic ram, and the entire process is continuously 
repeated until the error signal is reduced to zero. As the error signal 
decreases, the servovalve closes, control surface movement ceases, and 
equilibrium is established at the new position, 
a. Component Descriptions 

(1) Servovalve . The servovalve (Figure 2) used in a submarine 
steering control system is a two-stage, force feedback model, incorporating 
a torque motor actuated by an electrical error signal; a flapper-type 
pilot valve, the flapper itself also serving as the armature of the torque 
motor; and a spool valve, actuated by the fluid flow caused by the movement 
of the pilot flapper valve. Movement of the spool results in the direction 
of fluid flow to the actuating ram. 

(2) Actuator . The actuator is a linear, two-way hydraulic 
ram. Initially, it was assumed that the actuator was of compensated con- 
struction; that is, it was modeled as if it had been designed with equal 
effective areas on both sides of the piston. O-rings are incorporated 
into the design of the actuator, eliminating essentially all leakage 
from one side of the piston to the other. 
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FIGURE 2 
SERVOVALVE 



(3) Mechanical Linkage . The linkage connecting the hydraulic 



ram to the rudder consists of a long shaft, support bearings, a packing 
at the penetration of the pressure hull, a crosshead at the end of the 
shaft opposite the ram, and a rod connecting the crosshead to the rudder 
tiller. The crosshead provides a supported, pinned joint from which the 
connecting rod originates. It was assumed to be constrained from lateral 
movement by an open-ended sleeve which was flooded with sea water. 

(4) Rudder . The rudder is composed of two portions, an upper 
rudder and a lower rudder, interconnected to the tiller by a single shaft, 
and located on the top and bottom of the ship’s hull, respectively. The 
rudder is symmetric along the chord length, and is designed for operation 
within the range of left 35 degrees rudder to right 35 degrees. 

(5) Other Components . Other components include piping; an 
electric transmitter associated with the control wheel; an error feedback 
servo associated with the actuator; a hydraulic pressure supply which is 
assumed for design purposes to be constant in pressure, and not subject 

to influence by steering system conditions; and a petroleum-base hydraulic 
fluid. 

C. NSRDC BASELINE SIMULATION 

The following sections describe the work done at the NSRDC in order 
to simulate the submarine steering control system. The governing equations 
for each of the components are stated as they were derived by the NSRDC, 
along with any necessary assumptions. A schematic circuit diagram of 
the analog computer arrangement is presented, and comments from the 
NSRDC about the nature of the problem and suspected origins of the oscil- 
lation are stated. 
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1. Governing Equations and Assumptions 



The governing equations that were provided for each of the 
components are listed verbatim. Only the governing equations, assump- 
tions, and values for constants were communicated to the writer; the 
procedure by which the expressions were derived is unknown. See 
Appendix A for explanations of the notation used, 
a. Servo Valve 

The servo valve equations include expressions for the 
valve flapper displacement, spool valve displacement, nonlinearity com- 
pensations, and fluid flow rates through the valve. 

(1) Valve Flapper Displacement . The valve flapper 
displacement, x, is expressed as 



x =* sat 




(K 1 (0 O 



-6)-K 3 z)j 



+.003 



-.003 



( 1 ) 



where is the amplifier gain coefficient; is the flapper centering 
spring rate, 0^ is the desired position of the rudder; 6 is the actual 
position of the rudder, is the wand spring rage, and z is the present 
valve spool position. No explicit assumptions were given with regard to 
this equation. 

(2) Valve Spool Displacement . The valve spool displacement 
is obtained through the variation with time of the valve spool velocity z 



z * < 



0.5 



K.x 

4 



z >.5 and SGN(z) 
otherwise 



SGN(x) 



( 2 ) 



where is the pilot-stage valve flow coefficient. This expression is 
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best explained by stating that the spool velocity is proportional to the 
flapper valve displacement, until the spool is fully displaced against 
the valve body, where the movement must of course cease. The velocity 
again becomes non-zero only if the direction of movement is reversed, and 
the valve is allowed to travel away from the hard restraining block. 

Nonlinear effects in the valve, such as varying flow gain 
(i.e., flow is not proportional to spool displacement, due to the shape 
of the ports used in the spool valve) , spool friction, deadband, and torque 
motor hysteresis are included in the next equations. 

Dead band is included in the system through the implementa- 
tion of the following equations: 

f z - db where z >db, ^>0 

z =< z + db where |zl>db, z<0 (3) 

0 if -db<z<db 



where z* is the valve spool displacement after adjustment for spool dead 
zone. These expressions give the effective spool displacement, i.e., 
adjust the spool displacement to the value it would have for the given 
flow if no dead zone were constructed in the valve. 

Multiple flow gain is accommodated by combining z f with 

z *1 z 1 1 : 



z" = . 5z ' + . 5 z * I z'l . (4) 

where z" is the valve spool displacement after spool dead zone, and after 
multiple flow gain. 

Compensation for other nonlinearities such as spool 
friction and torque motor hysteresis may have been accomplished through 
adjustment of the coefficients in the equations for spool dead band, and 



1 / 



multiple flow gain. No additional expressions for these nonlinearities 
were given. 

(3) Flow Equation . The flow rate in the spool is related to 
the corrected spool displacement z n , the load pressure drop and the 
supply pressure P^ by the following expression: 



Q = z"K 5 sgn(P s -sgn(z)P L ) J P s - sgn(z)P L (5) 

where is the main-stage valve flow coefficient. No explicit assumptions 
were listed for this equation, 
b. Actuating Ram 

The position of the actuating ram is described in two ways; by 
the position of the rudder itself, and by the integral of the flow rate 
which has passed through the valve. 

(1) Ram Position in Terms of Rudder Position . The ram position 
y is related to the rudder position, with correction for mechanical resilience 
in the linkage, by 



y 



Lsin0 




( 6 ) 



where L is the length of the tiller, 0 is the rudder angle in degrees, 

is the mechanical resilience coefficient of the system, and A is the actuator 
piston area. 

(2) Ram Velocity Expressed in Terms of Pressure Change and Flow . 

The NSRDC determined that the ram velocity, flow rate, and load pressure 

have the following relationship to one another: 

Q = -Ay + k 6 p l (7) 

where K is the system hydraulic resiliance, and P, is the rate of change 
6 L 

of load pressure drop. 
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In the development of these equations it was assumed that 
there was no leakage from the control valve to the actuating cylinder at 
line connections, seals, and so forth, and that there was no internal 
leakage in the actuating cylinder itself. Also, it was assumed that 
the actuator used was compensated, or incorporated equal areas on opposite 
sides of the piston. 

c. Forces and Torques on the System 

All forces on the system were expressed in terms of the equiva- 
lent torques which were placed on the rudder structure, and a moment 
balance was written: 

16 » T + T + T_ (8) 

p n r 

where I is the polar moment of inertia of the rudder mass and the included 

water mass, 0 is the angular acceleration of the rudder, T is the torque 

P 

developed as a result of the pressure drop across the actuator, is the 

hydrodynamic torque, and T is the friction torque. 

r 

(1) Friction Torque . is the most complicated of the 

torques to describe* It is the combined friction torques from the shaft 
seals, the resultant of the normal force on the rudder bearing surface, 
and any other sources of sliding (coulomb) friction present in the system* 
Stationary objects in contact with each other stick together until sub- 
jected to forces larger than the friction force present between the 
objects. Additionally, the level of applied force necessary to sustain 
movement at very low speeds usually exceeds the force needed at higher 
velocities. In order to model these system characteristics, the NSRDC 

developed the following exptessions for T : 

r 
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r 



-< 



-C T „- T p3 if |i,*Tp|< T s 
-T s sgn(T H .T p ) if | T H *Tp) >Ts) 



6 < e 



( 9 ) 



V 



-T c sgn(0) 



lei 



>e 



where e is theoretically zero, but in practice is set to approximately 
0.2 percent of the maximum angular velocity of the rudder. = kTg, 
where k is slightly less than one, typically 0.95, and Tg is expressed by 



T s c d r 



IV + 






■ + K 1 1 0 +A 



Pi)' 



+ K 



10 



( 10 ) 



where is the coefficient of friction of the stock bearing, W is the 
weight of the rudder, R is the rudder stock radius, K 9 is the lift co- 
efficient of the rudder, is the drag coefficient of the rudder, and 

is the seal friction coefficient. It should be noted that the 

expression for T does not include any terms for damping proportional 
to rudder turning rate, 

(2) Hydrodynamic Torque . T is the torque generated by the 
rudder when it is placed at some angle of attack other than zero, due to 
the fact that the center of pressure on a rudder does not coincide with 
the turning axis of the structure. The hydrodynamic torque is expressed 
as 

v(^ (K 7 9 -V |el } (11) 



where V is the ship velocity, is the maximum ship velocity, and K 7 and 

K are hydrodynamic torque load coefficients. By adjusting the values of 
8 

and Kg, the hydrodynamic torque variation with rudder position can be 
adjusted to approximate almost any rudder design. If K., and Kg are 
redefined, the above expression can be written in nondimensional form: 
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HYDRODYNAMIC torque 




FIGURE 3 

HYDRODYNAMIC TORQUE BEHAVIOR 
FOR VARIOUS VALUES OF K y AND Kg 
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( 12 ) 



H 



1 Hmax 





0 _ 

0 

max 



where is the maximum hydrodynamic torque corresponding to Vq, and 

0 max is the maximum rudder angle. See Figure 3 for curves of T^/T Hmax 
for varying values of and K g . It can be seen that a variety of rudder 
torque characteristics can be simulated by adjustment of and Kg. 

(3) Pressure Torque . The pressure torque Tp is the torque 
that is placed on the stock of the rudder by the hydraulic ram, through 
the connecting linkage. It can be expressed as 



T p = LAPpCos 0 , (13) 

where L is the tiller length, A is the area of the actuator piston, and 
Pp is the load pressure drop. 

2. Analog Computer Simulation 
a. Hardware and Parameters 

Figure 4 is a schematic diagram of the analog computer program 
which was designed to model the governing equations. Table I gives a list 
of typical system dimensions and other parameters which were used to set 
simulation values. Table II gives a list of the potentiometer settings 
for the values listed in Table I, and also lists the expressions which 
were used to develop the settings. 

3. Discussion of the Baseline Simulation 
a. Performance Curves 

Examination of sample recordings obtained from the NSRDC 
showed that oscillations in velocity and pressure were present in the 
simulation outputs. The NSRDC indicated that the oscillation appeared 
to be a function of the inertia of the steering linkage, ram, and rudder, 
and of the system mechanical and hydraulic resiliences. Variation of 
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Table I 



Submarine Steering System Parameters Used in the Baseline Simulation 



Parameter 


Value 


Parameter 


Value 


A 


80 in 2 


L 


30 in 


C d 


.05 


p 

max 


6,000 psi 


db 


.05 


p 

Lmax 


3,000 psi 


I 


9100 in-lbf-sec 2 


P S 


3,000 psi 


K 1 


. 6 


Qmax 


278.86 in'Vsec 


K 2 


110 


R 


18 in 


K 3 


5.0 


T 

Fmax 


1 . 1661X10 6 in-lbf 


k 4 


375 


T 

Hmax 


3 . 93X10 6 in-lbf 


K 5 


7.20 


T 

Pm ax 


7 . 2X10 6 in-lbf 


K 6 


0.015 


W 


25,000 lbm 


K 7 


1.5 


X 

max 


.002 in. 


00 


2.5 


y 

y max 


17 in. 


K 9 


30,000 


z 

max 


.5 in. 


K 10 


150,000 


0 

max 


35 degrees 


K 11 


3,000 


0 

max 


10 deg/sec. 


K 12 


800,000 


6, time scaling factor 


1 
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I 



Table II 



Potentiometer Expressions and Settings for Baseline Simulation 





(See Fig. 4) 




Potentiometer 


No. 


Expression 


Setting 


1 


e o K l 

100K 2 x max 


.9546 


2 


e m K l 

100K-.X 

2 max 


.9546 


3 


Simulation 

Technique 


.1 


4 


^4 x max 

10z maxB 


.1500 


5 


K 3 z max 

lOOK^x 

2 max 


.1136 


6 


db 


.05 


7 


db 


.05 


8 


multiple flow gain 


.5 


9 


multiple flow gain 


.5 


10 


Qmax (1st footnote) 

10 w 


.406 


11 


AL ®max Ctt / 1 80 D C 2n d footnote) 
20 Kg, ^Lmax 8 


.4654 


12 


D 

r Lmax/p 

' r max 


.5000 


13 


^/Pmax 


.5 


14 


A ^Lmax 

K 12 y max 


.0177 


15 


L 


.8824 



2y 

J max 



25 



\ 




16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 



Expression 



Setting 



A^Lmax 

3.45X10 5 


(3rd 


footnote) 


K, ,0 
11 max 

3.45X10 5 


(3rd 


footnote) 


% 9 max 


(3rd 


footnote) 


1.075X106 






W 


(3rd 


footnote) 



1.075X10 6 



r 3.45X10 5 / A' 1 2 ( ‘ 3rd 

( '1.129X10 6 ) 

( 1.075XI06 /K129xl0 6 ) 2 C3a 



footnote) 



footnote) 



KlO 

T 

Fmax 

Cd R (1 . 129X10^) (3rd footnote) 
T 

Fmax 

e 



^Fmax 






^Frna y (T c/Ts) 

1000 IB 
max 

max 

®max 

®max (4th footnote) 

200 

^8 

10 



.6957 

.3044 

.9768 

.0233 

.6934 

.9066 
. 1286 

.8714 

.0010 
.1282 

.1269 

2857 

1750 

25 
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Potentiometer 

No. 



Setting 



Expression 

30 Ky >15 

10 

31 T Hmax .3370 

I® ^Fmax 



32 



Tpmax 

10 Tp max 



.6175 



1. This expression was later found to be incorrect. See the text, 

Initial Investigation. 

2. This expression was later found to be incorrect. See the text. 

Initial Investigation . 

3. Numbers appearing in the numerators or denominators of these expressions 
represent scaling values calculated by the NSRDC. 

4. Adjustment for SIN - COS generator design requirements. 
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these parameters modified the fundamental frequency of the oscillation, 
but did not eliminate it. It was found that the oscillation could be 
eliminated by the artificial introduction of viscous damping in the 
system, but only in much greater amounts than seemed to be present in 
the actual system, 

b. Suspected Origins of the Oscillations 

After some study of the baseline simulation, the NSRDC felt 
that the problem was located in either their recording equipment associated 
with their analog computer, in the mathematical models which had been 
developed, in settings of the simulator parameters, or in the basic assump- 
tions used in the derivations of the governing equations. 
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II . INITIAL INVESTIGATION 



The initial investigation was conducted in order to gain a more 
complete understanding of the simulation, and to determine problem areas 
requiring extensive investigation. The baseline simulation was duplicated 
on a large analog computer, techniques and models currently used in simulat- 
ing hydraulic control systems were studied, and information was gathered 
from personnel familiar with the baseline simulation development and sub- 
marine steering control systems. The baseline simulation was continually 
evaluated in terms of the information which had been gained, and several 
problem areas were found which required further analysis. 

A. DUPLICATION OF THE BASELINE SIMULATION 

The baseline simulation required duplication on local equipment in 
order to provide a check of the correctness of the programming completed 
at the NSRDC. The Naval Postgraduate School f s COMCOR CI-5000 analog 
computer was used in this effort. The CI-5000 is a general purpose 
computing system, including analog, digital control, and patchable logic 
sections, and incorporates enough components of various types to accommo- 
date the baseline simulation. 

The writer first attempted to directly adapt the baseline simulation 
analog schematic directly to the CI-5000, adjusting the simulation only 
to compensate for CI-5000 peculiarities, thereby assuming that the equations 
and potentiometer expressions of the baseline simulation were correct, and 
that the analog schematic diagram accurately modeled the governing equations. 
This first attempt was unsuccessful; analog components rapidly overloaded, 
and the system failed to reproduce the sample curves that the baseline 



29 



simulation at the NSRDC had produced when operated under identical 
conditions. Consequently, the baseline simulation’s governing equations, 
potentiometer expressions, and analog programming were checked for mathe- 
matical and functional correctness. 

The governing equations were found to be consistent with the sign 
conventions evidently adopted by the developer (See Figure 1), with the 
exception of equation (7), 

Q * -Ay + K & P l . 

Examination of the geometry of the system (Figure 1) revealed that the 
sign of the y - term should have been positive. Thus, equation (7) should 
nave read 

Q « Ay + . (7a) 

The NSRDC later confirmed that this was the case. 

Potentiometer expressions were found to be consistent with the 
baseline simulation schematic diagram and the governing equations, with 
the exception of potentiometers 10 and 11, which were associated with 
integrator 9 of the baseline simulation (Figure 4). A term in the denom- 
inator of each potentiometer expression apparently had been dropped; the 

2 

denominators should have read 10(P f )[K. + (A /K no )] instead of 10(K^P T ) 

Lmax 6 12 6 Lmax 

and 20(K.P f ) for pots 10 and 11, respectively. The need for the 

o Lmax 

additional term in each expression’s denominator was later confirmed by 
the NSRDC. The additional terms modified the vaLues of the two potenti- 
ometers significantly; the value of pot 10 changed from 0.6197 to 0.4060; 
the value of pot 11 changed from 0.4654 with a gain of 20 on integrator 
9 to 0.6060 with a gain of 10. Tests later showed that the omission of 
these coefficient terms was the primary reason for the initial failure 
to duplicate the baseline simulation. 
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Except for the two potentiometer settings, the analog program (Fig. 4) 
appeared to accurately represent the verified governing equations. Accord- 
ingly, another attempt was made to adapt the program to the CI-5000. Only 
slight deviations from the baseline schematic diagram were necessary in 
order to accommodate the CI-5000 component peculiarities. 

When the adapted baseline simulation was operated under the potenti- 
ometer conditions listed in Table II, with the exception of pots 9 and 10 
which were set as previously modified, the system behaved in the same 
manner as it had during operation at the NSRDC. Oscillations identical 

to baseline forms were observed in 0 and P (Figure 5). As a result of 

L 

this duplication of performance, using several different types of recording 
equipment that exhibited varying response characteristics, it appeared that 
the recording equipment at the NSRDC was not the source of the oscillations. 
Further, since all changes in mathematics and in the settings of the simu- 
lator parameters were implemented before the baseline performance was 
successfully duplicated, it appeared that the errors which were found were 
generated through the communication of the baseline simulation to the 
writer, and were not present in the program used to generate the sample 
curves. Apparently, then, the mathematical models were consistent in 
sign conventions and had been manipulated correctly in the development of 
the analog board, the potentiometers were correctly scaled, and the analog 
wiring was correct. Thus, the simulation appeared to be correct within 
the framework of the fundamental assumptions involved in the mathematical 
formulations; however, the oscillations were still present in the output. 

Two situations seemed to be possible. First, perhaps the oscillations 
evidenced in the NSRDC data were not representative of the behavior present 
in real submarine steering systems, and the simulation was therefore 
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FIGURE 5 

BASELINE SIMULATION RESPONSE 
e o = 35 DEGREES 
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inadequate. This situation could have occurred because omissions were 
made during the formulation of the baseline simulation governing equations , 
such as a miscalculation of the significance of some physical process 
present in a steering system. The situation could also have occurred 
because system materials, dimensions, and properties were incorrectly 
represented in the development of the numerical values using in the "typical" 
system simulation, due to the inclusion of inconsistent units conversions, 
unrealistic dimensions, or faulty physical materials constants. 

Second, perhaps the oscillations appearing in the NSRDC data were 
representative of real system behavior, and the difficulty had its origins 
in the harboring of an incorrect knowledge of the actual performance of a 
ship's steering control system. No actual system recordings were submitted 
with the baseline simulation, and therefore no immediate comparison was 
possible between the simulation performance and real conditions in the 
system. 

B. BACKGROUND INVESTIGATION 

1. Study of Currently Used Models 

Assuming that the oscillations observed in the baseline simulation 
were not representative of real steering system behavior led to a back- 
ground investigation of techniques and models currently used in simulating 
hydraulic systems in order to develop a clear knowledge of the fundamental 
assumptions made during the development of the baseline simulation. This 
background examination consisted of studies in the hydraulic controls 
field, and communication with personnel familiar with the simulation and 
ship steering control systems. 

The search for currently used models and techniques produced several 
results. It was noted that if the expressions for the baseline simulation 
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were transformed into Laplace notation, and deadband and multiple flow 
gain effects were temporarily neglected, the resultant expression would 
be 

V K 2 

Z(s) = = — E(s) (14) 

S+K 3 /K 2 

where Z and E represent the Laplace-domain variables for z, the valve 
spool position, and e, the error signal, respectively, and s is the 
Laplace transform variable. This function was identical in form to one 
developed by Morse [lj for a valve incorporating centering springs instead 
of force feedback, for use in systems where the load was assumed to bu 
very small compared with the available hydraulic power. The actual valve 
being used, however, consisted of a two-stage, force-feedback design being 
used to control a load of the same order of magnitude as the available 
hydraulic power. Merritt [2] developed a transfer function block diagram 
for a two-stage, force-feedback servovalve, which incorporated the general 
form shown in Figure 6. A vast degree of simplification of Merritt's 
expression obviously would have been necessary in order to arrive at the 
expression used in the baseline simulation. Quite possibly, the coef- 
ficients indicated by Merritt in many instances were of such magnitudes 
that their effects were not important in the actual performance of sub- 
marine steering system servovalves, since high frequency effects would 
become unimportant when other system frequencies were quite low. Further- 
more, the pressure feedback to the flapper valve that Merritt included in 
his expression was possibly unimportant in this particular case. However, 
it seemed advisable to check the numerical value of each term in the 
detailed expression, in order to evaluate the validity of the baseline 
simulation's expression. 
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FIGURE 6 

BLOCK DIAGRAM OF A TWO -STAGE SERVOVALVE 

WITH FORCE FEEDBACK 




It was found [3] that many of the servovalve design values, needed to 
calculate the coefficients in the Merritt block diagram, were difficult 
to obtain. It was also found that the expression used in the baseline 
simulation was considered to be sufficiently accurate for the particular 
design studies in which it would be used. Finally, it was determined 
through observations of the performance of the adapted simulation program 
that the oscillations did not appear to be affected by the servovalve 
action. Thus, further investigation of the validity of the servovalve 
simulation, did not appear to be justified until the origin of the oscilla- 
tions had been fully investigated. 




Figure 7. Functional Block Diagram of a Submarine Steering System 

Figure 7 is a function block diagram of a submarine steering 
control system. The control system was not entirely closed loop in con- 
struction, since the error signal controlling the function of the servo- 
valve was generated by the ram position rather than by the actual rudder 
position. This placement of the feedback element was consistent with 
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FIGURE 8 

CURVES OF SEVERAL PARAMETERS 
OF Cl -5000 SIMULATION 
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FIGURE 9 

“FRICTION LES$>' LOAD CONDI TION,T F = 0. 
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established design practice (4), since the inclusion of the compliant 
structure of the rudder and linkage within the servo loop could result in 
significant instability problems. Additionally, an electrical feedback 
element located outside the pressure hull of a submarine, in a location 
flooded with sea water, would pose several design problems, including the 
need for increased reliability due to its inaccessability when at sea. 




Figure 10. Circuit Modification of Fig. 4, 

Introduction of Viscous Damping 

p . , 0 6 P L T H z" y . z 

Examination of the curves for — , , and — 

o 6 r T z y z 

m m Lm Hm ra *ra m 

(Figure 8) generated by the CI-5000 simulation revealed that the oscilla- 
tions were present only in the linkage and the load portions of the simu- 
lator. Manipulations of the baseline simulation friction coefficients 
failed to eliminate or modify the effect, although total exclusion of any 
friction in the system resulted in only very small undamped oscillations 
in 9, and noticeable oscillations in 0 and P after the control surface 

L 

had essentially reached the desired position (see Figure 9). Further 
tests showed that small values of negative feedback across integrator 25, 
which were introduced through the use of the circuit modification indicated 
in Figure 10, resulted in elimination of the oscillations. This feedback 
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across the integrator would be representative of damping in an actual 
system containing such elements as dashpots. With no other friction in 

t 

the system, a feedback level of 0,005 - — was sufficient to eliminate the 

h 

* m 9 

limit cycle. Furthermore, feedback levels of 0,12 t— were sufficient to 

^m 

approach a critically damped condition in the load and linkage simulator, 

(See Figures 11 and 12.) If friction levels present in the baseline 

simulation were programmed into the CI-5000, approximately the same value 

of h - damping as in the frictionless configuration was required to produce 

critically damped behavior, while the sustained limit cycle was eventually 

eliminated by the presence of the baseline friction itself, 

# 

Damping terms proportional to & were totally absent from the baseline 
simulation equations. Intuitively, however, some damping would be present 
in any system in which viscous fluids such as sea water are displaced by 
moving objects, since work must be expended to keep objects moving at a 
constant rate, whether linearly or rotationally , through a real fluid. 

Thus, the neglect of damping in the linkage and rudder portion of the 
simulation did not appear to be realistic. A superficial examination of 
the system did not reveal any large sources of damping, or even sources 
which would provide a fraction of that needed to approach a critically 
damped condition in the linkage. Anticipating that an investigation of 
system damping could become quite involved, the writer identified this 
area as requiring first priority for further investigation, and postponed 
further analysis of it. 

Examination of steering control system drawings for the SSN-671 class 
submarine revealed that functional concepts of several features of the 
design had been altered in the process of the development of the idealized 
model. A compensated actuator had been assumed in the simulation devel- 
opment; i.e., equal areas had been assumed on both sides of the actuator 
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FIGURE I I 

VISCOUS DAMPING LEVEL . 005#-, 

Tp = 0. °max 
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FIGURE 12 

VISCOUS DAMPING LEVEL .I2§r ,T F =0. 

°max h 



42 



ram (Figure 1), while examination of the drawings revealed that the 
actuator on a typical submarine was usually uncompensated. It was 
learned that the effect of neglect of this departure from actual system 
design was thought to be small, and not important in terms of the projected 
simulation use [5], 

Further study of the drawings revealed that the crosshead configu- 
ration had been simplified. The crosshead, originally assumed to operate 
in a sleeve open at both ends, Figure 1, was in reality supported by a 
sleeve open at only one end in most submarine designs, Figure 13. Further 
investigation of the crosshead influence on the system seemed to be justi- 
fied, since the assembly resembled a dashpot in appearance, and possibly 
was capable of providing the source of damping which seemed to be needed 
in the system. Consequently, a detailed analysis of the effect of a 
typical crosshead structure on steering system performance seemed to be 
in order, and was identified as an area for further work. 

2. Verification of Typical System Values 

Several nondimensional coefficients, load parameters, and properties 
had been used to develop the simulation, such as the hydraulic resilience 
coefficient, the rudder mass, and inherent materials properties. Inaccurate 
development of system numerical values could have contributed to the 
system oscillations. Attempts were made to determine the values of many 
of these constants from basic principles, but the effort was thwarted 
due to a lack of information regarding dimensions, physical constants, 
system construction details, and algebraic expressions for any of the non- 
dimensional constants through This possible source of inaccurate 

system simulation was identified for further investigation, and was set 
aside. 
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3. Actual System Performance 

Although actual system recordings of 0 , 6 , and P L were mentioned 
[6] as not demonstrating the oscillations found in the simulation, such 
recordings were not forwarded to the writer. Because an accurate simula- 
tion of a system could not be claimed until it had been favorably compared 
with the actual system performance, it was felt that such a comparison 
would be beneficial, and the subject was identified as requiring further 
investigation. 

C. SUMMARY OF THE INITIAL INVESTIGATION 
1. Conclusions 

As a result of the duplication of the baseline simulation, and 
subsequent experiments performed with the program, the following conclu- 
sions were reached: 

(1) The recording equipment at the NSRDC had accurately recorded 
the performance of the simulation, and was not the source of the oscil- 
lations . 

(2) The equations and potentiometer settings were consistent 
within the framework of the fundamental assumptions involved in the 
simulation's mathematical development. 

(3) The oscillations were confined to the load portion of the 
simulator, and did not appear to be sustained by overall system instability. 

0 

(4) A negative feedback level of 0.005 t— was sufficient to 

. *m 
0 

eliminate the limit cycle oscillation, and 0.12 — achieved a critically 

0 

damped condition, in the load simulation under otherwise frictionless 

0 

conditions. A negative feedback level of 0.12 - — was also sufficient to 

9 

ra 

create a critically damped condition in the load simulator under baseline 
friction conditions. 
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Studies, and communications with personnel familiar with the 
simulation and with actual submarine steering control systems, led to 
the following conclusions: 

(1) The servovalve equation used in the baseline simulation was 
not designed to completely model either the type of servovalve used in 
the real system under study, or to be used with the type of load present 
in the actual system. The net effect of the substitution of a simple 
expression for a complicated one was not immediately determined, except 
that it did not appear to be capable of affecting the observed simulator 
load oscillations. 

(2) The actuator had been modeled in the system simulation as 

a compensated ram, rather than being modeled as the uncompensated actuator 
used in real systems. Although the effect of this simplification on the 
system was unknown, communication with experienced personnel seemed to 
indicate that the effect was unimportant. 

(3) The actual crosshead design had been modeled in such a way 
that damping effects possibly present in the component were not taken into 
account* Further study seemed warranted. 

(4) Viscous damping effects in the load had been neglected in 
the development of the simulation, although such effects controlled the 
unwanted system oscillations, and normally would be present to some degree 
in any system incorporating the movement of structures through real fluids. 

A search for viscous damping sources in the load was assigned first priority 
for further investigation. 

(5) The origins of many of the nondimensional constants used in 
the simulation were not obvious, and efforts to verify their correctness 
were hindered by the lack of their related algebraic expressions. 
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(6) The realism of the simulator performance needed to be 
checked through the comparison of simulator and actual system recordings. 
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HI. LOAD AND LINKAGE DAMPING 



Of the areas identified in the initial investigation as requiring 
further research, only the search for sources of linkage and load damping 
was pursued in the present investigation. Analysis of the submarine 
system identified three locations where damping might possibly be present: 
the system’s dynamic seals and bearings, the crosshead structure, or the 
rudder itself. No other components appeared to be designed in such a 
manner that they could possibly produce an appreciable quantity of damping. 

In order to determine levels of damping in the various components, 
detailed information was needed concerning dimensions, clearances, materials, 
constructions, and lubrication for the linkage parts and rudder. The writer 
was referred [ 7 ] to Mr. Leeland Smith [ 8 ] of the Mare Island Naval Shipyard, 
who provided detailed drawings of the steering control system for the SSN 
671 class submarine. 

A. BEARINGS AND DYNAMIC SEALS 

The possibility was considered that the expressions employed for the 
friction present in the bearings and dynamic seals (equations 9 and 10) 
were not representative of the actual processes present in the submarine 
system, and that the inherent friction of the bearings and seals was 
better described as displaying viscous behavior. 

Keller [9 ] indicated that compression packings, such as those used 
at the pressure hull penetration and the actuator-piston rod interface, 
are normally lubricated by allowing very slight amounts of leakage past the 
seal. Examination of the SSN-671 drawings revealed that the pressure hull 
packing was additionally lubricated with grease. Despite these indications 
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that viscous effects could be present in these seals, discussions with 
Mr. Smith, William Kendall [10], and Leonard Chase in] of the Mare Island 
Naval Shipyard resulted in the writer's learning that friction in such 
seals was best approximated as sliding friction, since friction loads were 
normally observed to be fairly independent of actuator rod velocity, due 
to the velocities and allowable leakage rates present in the system. Thus, 
the baseline model of seal friction appeared to be justified. 

Study of the rudder bearings seemed to support the use of the baseline 
simulation's bearing friction expressions. The large rudder weight, com- 
bined with hydrodynamic loads, small rates of rotation present in the 
rudder stock, and the ill-defined lubricating qualities of sea water, 
resulted in enough analytical uncertainties being present to discourage 
any attempt to model such bearings as either slipper-type or journal- type 
surfaces. Again, the baseline simulation expression seemed to provide a 
bearing surface description that probably was as accurate in magnitude and 
parameter dependence as any other expression that could be analytically 
developed using known conditions and existing lubrication theory. Experi- 
mental and analytical work of a scale beyond the scope of 'this project 
could prove useful in future design work involving bearings and seals of 
the types present in the submarine system. 

B. CROSSHEAD 

1 . General Description 

Examination of the SSN-671 drawings, and discussions with Mssrs. 
Smith, Kendall, and Chase concerning general construction details of 
crossheads, revealed that, as a rule, crossheads somewhat resembled dash- 
pots in construction, due to their piston-in sleeve arrangements, and their 
free-flooded nature. Provisions were normally made in the construction of 
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crossheads to eliminate any possible damping effect by the drilling of 
large holes through the piston portion of the structure. Because damping 
sources were being sought, the writer decided to analyze the crosshead in 
order to determine what damping effects, if any, were present. 

Figure 13 shows a crosshead design based on information taken from 
the SSN-671 drawings. The crosshead consists of a partially bored out 
cylinder 19 inches in length, with a diameter of 11.986 inches, and four 
one-inch I.D., eight-inch in length holes axially bored through the solid 
head of the piston. This piston rides in a twelve-inch I.D. tube, typically 
sealed along its length except at the rear end, which is open to allow the 
passage and free rotation of the link connecting the actuator shaft ex- 
tension and the rudder tiller. The piston incorporates several lock- 
preventing 0.75 in. X 0.125 in. circumferential grooves, and two 0.125 in. 

X 0.5 in. axial slots along the outside radius. As a result of the radial 
clearance between the piston and the sleeve, the slots, and the holes in 
the crosshead, the entire sleeve is free flooded, and flow of sea water 
past the crosshead must take place as it is moved in and out of the sleeve. 

2 . Assumptions 

In approximating the behavior of the flow past the crosshead and 
the viscous drag associated with that flow, it was assumed that the fluid 
velocities through the holes, slots, and radial clearance of the crosshead 
were considerably greater than the motion of the piston relative to the 
sleeve, thereby implying the existence of fixed orifices of irregular 
shape. Laminar flow was also assumed. 

3. Equations 

a. Eccentric Piston in a Sleeve 
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For an annulus formed by a solid eccentric cylinder in a sleeve 



with no relative movement involved, the flow rate can be expressed as [12] 



Q c = 



iTrc 

6 ITT 



m 



( P l- P 2) 



(15) 



where r is the radius of the sleeve, y is the absolute viscosity of the 
fluid, 1 is the length of the tube, e is the eccentricity of the piston with 
respect to the center of the sleeve, and (P]_-P2) is the pressure drop across 
the length of the tube. 

b. Laminar Flow Through Rectangular Passages 
For the slots [13] , 

3 



Q, = 



wh‘ 

12yl 



, 192h . , irw 

1 " 1— tanh 2h 

7T W 






- P 2' 



(16) 



where w is the width of the slot, h is the depth, and all other symbols 
retain their previous meanings. 

c. Flow Through Short-Tube Orifices 

Merritt [13] • combining the work of Langhaar [ 14 ] , Kreith [ 15 ] , 
and Shapiro [j^] , developed an expression for the flow rate through short 
tube orifices : 




(17) 



(18) 



where C is the loss coefficient, p is the fluid density, D is the hole 



diameter. Re is the Reynolds Number based on the mean velocity of the 
flow in the tube and the hole diameter, and all other variables retain 
their previous meanings. If equations 15, 16, and 17 are manipulated 
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briefly, the pressure drop Pj-P 9 can be expressed as a function of the flow 
rate Q q , the kinematic viscosity v, D, and p: 



P l- P 2 




1.5+6.87j™I 0^12.5^1 

2.28+ 167rvl C^<12.5™1 
Q 0 



(19) 



If a total mass balance is made across the crosshead, it is found that 



Aa-i y = N o Q o + N s Q s + Q c (20) 

where Aqj is. the cross-sectional area of the crosshead, n o is the number of 
holes in the crosshead, and N s is the number of slots. 

The force caused by the pressure difference across the crosshead 
can be expressed as an equivalent torque at the rudder pivot, by multi- 
plying the crosshead force F = by the length L of the tiller 

arm : 

T CH = A c„ 1 ( p r p 2) (21) 

Figure 14 shows curves of Tqj versus y generated through the solution of 
equations (18) ~(20) using the previously listed crosshead dimensions, p 
equal to 64 lbm/ft^, and v equal to 1.41 X 10~^ ft 2 /sec (water at 50 degrees 
F) , for various sizes of piston holes, assuming four holes and two slots, 
and an eccentric piston in the sleeve. To generate the curves, equations 
(18) and (19) were solved for P^-P2 as a function of Q^. Then and 

were calculated, and appropriate multiples of Qq, Qg, and Q^, were used in 
accordance with equation (20) to produce y. Equation (21) was then used to 
determine T^ 

For any crosshead incorporating four holes any larger than 0.1 
inches in diameter, it was found that zero pressure differential existed. 

As the diameter of the holes was decreased, with the dimensions and number 
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FIGURE 14 

CROSSHEAD TORQUE 
FOR VARIOUS HOLE SIZES 
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of slots and the parameters for the annulus remaining constant, a pressure 
differential began to exist across the structure. Finally, as the holes 
were eliminated, approached levels sufficient to create significant 

values of Tq^. For the values of the zero-holes curve used in Figure 14, 

T ch c ould be expressed as 



T Ch = 7.48X10 3 y 



y e 

or, since - — ro -r- 

0 



7 m 



max 



t CH = (7 .48X10^ymax) (5 



max 



( 22 ) 



(23) 



If the above equations were scaled for introduction into the 
simulation, potentiometer 33 would be set according to the expression 



7.48XlCrymax 

= t 

33 100 I 0 max 

For I = 9100, and y max = 5 inches per second. 



(24) 



a 33 = .00412. 

Since this value of the potentiometer represented the maximum available 
T ch , it was seen that the actual effect of the crosshead on the system 
damping was very small. In fact, using one-inch diameter holes would make 
the crosshead damping essentially zero. However, since a feedback level of 

0 

.005 *r — was previously found to be sufficient to eliminate limit cycles in 
e m 

0 

the "frictionless" system simulation, .00412 _ — could assist other small 

6 m 

sources in the damping of the oscillations in the baseline simulation. Al- 
though the results of the investigation of the crosshead damping were less 
spectacular than had been desired originally, since the 0- feedback sufficient 
to create a critically damped condition in the load had not been justified, 
the author felt that the identification of a small source of damping had 
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been achieved, that identification possibly proving to be of some worth in 
future design work if a small amount of damping somehow proved necessary. 



could create more damping. However, for this to be done, the system would 
need to be sealed at both ends, due to the approach of P ? to the vapor 
pressure of water in the contained volume when the crosshead was moving 
out of the tube. Of course, the open end of the crosshead sleeve was not 
capable of being closed, due to the primary function of the crosshead, 
which was to provide a portion of a slider-crank mechanism. The previously 
developed levels of Tqj thus appeared to be the maximum practicable. 
Additional damping could be generated, if needed, through the use of a 
large dashpot in parallel with the actuating linkage. However, such a 
dashpot was not apparently needed, as was shown through investigation of 
the inherent damping of the rudder. 

C . RUDDER 

The description of the moments involved in the rotation of a rudder 
about its pivot, while moving through a fluid, would in general include 
functions involving the zero, first, and second-order derivatives of rudder 
angle with respect to time. That is, if a moment balance were written for 
a rotating rudder, the equation could in general take the following form: 



Further consideration of the natures of f^ and f revealed that expressions 
for these two functions were present in the baseline simulation. If f^ = 



Some argument could have been made that further modification 



m = f : (0) + f 2 (e) + f 3 (e) 



( 25 ) 





( 26 ) 



q 5 



( 27 ) 



and f is recognized as being of the form 

h ■ » 

where I is the polar moment of inertia of the rudder and its fluid added 
mass. If f 2 = f 2 = 0 , fj is recognized as the sum of all 0 - dependent 
parameters in the rudder, including the hydrodynamic torque caused by the 
rudder lift and the difference between the location of the rudder’s pivot 
and its center of pressure, and the 0 - dependent portion of the system 
friction. 

The function f o (0) was assumed equal to zero in the baseline simulation; 
that is , no e - dependent effects were included in the rudder’s equation of 
motion (Eqn. 8). The validity of this assumption was examined in the final 
portion of the investigation of load and linkage damping. 

1. Existence of a Non-Zero Function f 2 ) 

Initial considerations of the problem involved the search for any 
indication that a non-zero function f 9 (0) even existed. Intuitively, such 
a relation should exist in some magnitude, at least under a condition where 
the rudder was subjected to a constant turning rate while immersed in a 
stagnant, viscous fluid, since power would necessarily be expended to keep 
such a structure moving under these conditions, thus resulting in some level 
of torque being applied to the rudder shaft. This applied torque, then, 
would by Eqn. (25) be equal to f 9 (0), and therefore f 2 could not be identically 
zero. 

Another indication of the existence of a non- zero f2 was provided 
by Landweber [17], who indicated that forces on lifting surfaces in unsteady 
motion were due to a combination of inertia effects and circulation, and 
that a transient variation did exist in the circulation around a thin air- 
foil impulsively given an angle of attack a, causing a transient behavior 
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of the lift coefficient of such a surface. Additionally, Theodorsen [18] 
studied the forces on a foil undergoing harmonic oscillations in connection 
with studies of control surface flutter and aerodynamic instability, and 
found that the actual lift coefficient was related to a quasi-steady lift 
coefficient corresponding to the instantaneous angle of attack of the foil, 
and time. Thus, the existence of a non-zero function f 2 ^ O ^ ^ w ^il e not 
confirmed by the above references, was at least shown to be plausible. The 
writer therefore attempted to locate any work that had been done in the 
analysis of the forces and turning moments on airfoils or rudders subjected 
to constant turning rate. 

2 . Damping in a Rudder Immersed in a Stationary Fluid 



by a rudder subjected to a constant rate of rotation while immersed in a 
viscous fluid. Accordingly, the writer attempted to approximate the drag 
which could be expected, using elementary concepts. 



No references were found describing the resisting torque developed 




dr 



\ 



s 




Figure 15. 



Rudder Construction Assumed For 
Stationary Rudder Torque Calculation 
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The rudder was approximated by a series of very small flat plates, 
assumed to be separated from each other so that the flow could move around 
the plates (Figure 15). For a flat plate placed in a stream flowing perpen- 
dicular to the plane of the plate, the drag force F^ can be expressed as 



F 



D “ 




A 



P 



( 28 ) 



where in this case is the drag coefficient, p is the density, is the 
area of the small plate, and the velocity V is related to the radius r and 
the angular rate of rotation 0 by V = r0. Thus, for a differential area 
dA n , 



dF D ' C d 






dA. 



If the area is assumed to be equal to the span of the rudder S times the 
incremental distance dr,dAp = Sdr , then 

(4 i) 



• 2 2 

dF D = C D [ P V - I Sdr 



The incremental torque dT^ caused by F^ developed at a distance r from the 
pivot point of the rudder would be 

: 2 



dT R “ C D ( T" ) Sr3dr 



(29) 



Thus, for the portion of the plate from r - 0 to r - r Q , 



W -2 

t r = — 1— 8 



( 30 ) 



For a flat plate in a normal stream with a large S/dr ratio, C D is 
approximately 1.95 [ 19 ] - For the two rudders of a typical submarine, 
the following physical parameters are representative: 
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Parameter 



Value 



s 




10 


Feet 


V 


rear portion 


7 


Feet 


r o> 


forward portion 


3 


Feet 


p. 


sea water 


64 


lbm/ft 



The entire 0 - dependent torque would be equal to the sum of the effects 
of both upper and lower rudders, including forward and aft portions: 

T 0 ” 2 ( Tr forward + Tr aft^ • ( 31 ) 

Substituting in the previously mentioned values, it is found that 

Tq = 88 0 2 , in. -lb f (32) 

when 6 is expressed in degrees per second. If e = 10 deg/sec, T- = 8800 

emax 

in. -lbf. If an asymmetric square analog component and a potentiometer were 
placed in the simulator as in Figure 16, this source of damping could be 




Figure 16. Analog Circuit Diagram for Introduction of T 0 

9 

The appropriate potentiometer setting expression for simulating T 0 would be 



a 34 



T 9 m 

1000 18 
m 



(33) 



59 






Substituting numerical values, 

a 34 = .968 X 10" 3 . (34) 

Thus, the initial estimate of available rudder drag at zero velocity 
indicated negligible amounts of damping were present in the rudder 
turned at a constant rate in a stationary fluid. The writer felt that 
the results were inconclusive, however, since the assumptions which were 
made in the calculation of T^ were very rough indeed. The rudder was 
not formed of many small strips with spaces between them; nor was the 
flow as simple as that described by the model. Such a complicated flow 
as would be present along the surface of a large flat plate rotating in 
a quiescent fluid required a depth of analysis beyond the scope of this 
study. 

3. Damping in a Rudder Immersed in a Moving Stream 

A literature search revealed that an analysis had been made of the 
torque behavior of a turning rudder immersed in a moving stream. Okada ( 20 j 
published a paper dealing with the effect, which was subsequently expanded 
by Kennard and Leibowitz [21] in their studies of rudder performance as 
affected by ship turning, steady change of rudder angle, and propellor race. 

Kennard and Leibowitz [22] noted that the twisting moment observed on 
a rudder shaft while the rudder is being turned may differ considerably 
from its value when the rudder is stationary. They went on to explain that 
this effect was due to circulatory motions around the foil caused by 
viscosity effects of the water, which subsequently generate vortex sheets 
in the water downstream from the rudder’s trailing edge, and cause moment- 
modifying variations in the pressure distribution over the rudder’s surface. 
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Variations from steady state twisting moments by as much as thirty percent 
were mentioned as being observed in rudders turning at constant rates. 

The mathematical difficulties associated with the problem were 
simplified by several approximations. The rudder was replaced by the flat 
plate, shown in Figure 17, of uniform chord C and infinite vertical span. 
The water was assumed to approach the plate at a uniform rate and at an 
angle of attack a; its velocity was denoted by U. The rudder twisting 
moment per unit span was identified as M^, and was assumed positive when 
operating in a direction to decrease a. The pivot point of the plate was 
identified as H, and was located at a distance hC from the leading edge, 
h being assigned some numerical value between 0.25 and 0.50. The rudder 
turning rate was assumed positive for increasing angles of a. 




Figure 17. Flat Plate Approximation to a Rudder 
Turning at Constant Rate 

M^, the twisting moment per unit span, was placed in nondimensional form 
by dividing it by (1/2) C 2 U 2 , to form the dimensionless moment coefficient 
C An expression for was then derived, with being represented as 
the sum of a quasi-stationary part Cj^q, and a 0 - dependent part, AC MH - 
Thus , 

C MH = C MH0 + AC MH (35 ' ) 
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FIGURE 18 
A VERSUS fm 
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where 



and 



C MHO = (4 _h ) 2 ttoi 



(36) 



A c mh 



( '-A) 




-h) ( 1 — 2h ) - (d. -h)A 



‘1 



(37) 



_ C da 
^ U dt 



and A is a function of XI and a . The effect of A is to increase the value 

of AC^jh as a increases, as is shown in Figure 18. 

Kennard and Leibowitz reported that the expression for C was not 

MHO 

entirely correct, due to the neglect of a nonlinear effect in its develop- 
ment. That nonlinear effect concerns the aft movement of the center of 
pressure in a rudder as a is increased, which causes a change in the sign 
of increased a. Also, they indicated that values for AC^ obtained 

from model tests suggested that the calculated values of AC^ were somewhat 
too high, although the scatter of the experimental data in the tests was 
considerable. They further concluded that the general trends resulting 
from turning real rudders in water were modeled in their theory, although 
values for AC^ were somewhat too high, such as a calculated increase in 
one case where h = 0.3 of roughly 48 percent, where in practice the increase 
due to AC^j ranged from 15 to 40 percent. It was suggested that tests with 
larger rudders would be worthwhile, and that refinement of the analysis was 
needed, expecially in including the effects of finite rudder spans. 

4 . Application to the Simulation Problem 

The analysis just described seemed to be applicable to the simulation 
problem. Although the expression for C^q was apparently unrealistic, the 

development ofAC was apparently useful in predicting the effect of rudder 
MH 

turning on the observed moments on a rudder, although calculated values 
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seemed to be somewhat higher than values which were experimentally obtained. 
In spite of its inherent shortcomings, incorporation of this modeling of 
the effect of rudder turning rate on the observed twisting moment seemed 
justified in light of the reported magnitudes of change in twisting moment 
observed in practice. 

In order to apply the above analysis to the simulation, several 
assumptions were made. It was assumed thatawas identical to 6 ; that is, 
that the fluid flowed parallel to centerline of the submarine. Analysis of 
drawings of the SSN-671 submarine rudder revealed that h was approximately 
equal to 0.28; in order to compensate for the reported high nature of the 
calculated AC^, and to simplify calculations, h was assumed to be 1/4, 
thus eliminating dependence on A. Also, the sign convention of was 
reversed to conform to the simulation^ convention. The above assumptions 
when used in the expression for AC^ resulted in 

(38) 

Re-dimensionalizing to form AM h , the change in twisting moment per unit 
span due to 0, and then multiplying by the span length S, resulted in the 
expression 

AM h 5 = - pc 3 sve . (39) 

For a SSN-671 rudder, the average value of C is about 10 feet, and S is 
also about 10 feet. In a typical submarine V can vary from zero to in 
excess of 35 ft/sec, and 0 was assumed in the baseline simulation to vary 
from zero to 10 degrees/sec. Assuming p = 64 lbm/ft 3 , the maximum value of 
the twisting moment, (AM^S) max based on Vo = 35 ft/sec, would be 

= 5.72 X 10 5 in-lbf . 

max 



(AM h S) 



(40) 



Typically, a submarine has two rudders, of approximately the same design. 
The total twisting moment available would therefore be twice the above 
value, or 1,144 X 10 in-lbf. An expression for AM^S/ (AM^S) would be 



am h s 



(AMyS) 



max 



■iyiu 



(41) 



AMrS/ ( is dependent only on for a given value of V, and therefore 

could be simulated as damping across integrator 25 as in Figure 10. The 
appropriate expression for the potentiometer setting a^ would be 



33 



( M S) (V/V ) 
H max o 

10 (100)B 

max 



(42) 



For V - V^, and the conditions otherwise assumed above, since V appeared to 

be set equal to for the other V - dependent parameters in the baseline 

♦ 

0 

simulation, a^~ « 0.1259. A level of 0.12 — feedback was previously found 
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to be sufficient to approach a critically damped behavior of the simulation, 
and therefore eliminate the oscillations from the system. Apparently, then, 
a theoretical justification for the inclusion of that much damping in the 
system had been found. The fact that AC was dependent on V did not allow 
the extension of the above damping levels to lower velocity, since as V 
decreases, so must the level of damping. The authors felt that the damping 
effect somehow assumed a significant value even at low or zero velocities, 
thus indicating the existence of some unexamined or undiscovered contributing 
effects such as those previously discussed, which would be independent of V. 
More work in the area of rudder twisting moment and lift variance with rudder 
turning angle and other parameters obviously would be useful. 

5. Summary 

The investigation of linkage and load damping resulted in the 
location of a significant source of 6 - damping in the motion of a rudder 
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immersed in a moving fluid. A significant amount of damping was evidenced 
by the large variation of the observed twisting moment produced by the 
turning of a rudder from that which would be observed if the rudder were 
held stationary at a given angle of attack. On the basis of rudimentary 
analysis, no large amounts of viscous damping were located in the packings, 
or the crosshead structure, although a slight amount of damping could be 
caused in the crosshead through redesign of the structure. A very rough 
estimate of the twisting moment observed on a rudder immersed in a stationary 
fluid indicated that little available damping existed unless the rudder is 
moved through the fluid. More analytical and experimental research needs 
to be done concerning the transient lift and twisting moment characteristics 
of moving and stationary rudders in fluid flow fields. 
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IV. RECOMMENDATIONS FOR FURTHER WORK 



Several topics were identified as requiring investigation during the 
progress of this project. Due to time limitations, none of them were 
pursued in detail. Subjects uncovered in the initial investigation which 
remained significant at the end of the period that this report covers 
included the need for a comparison of the actual and simulator performances 
under varying conditions, including a wide range of ship velocities. 

Studies of. the effect of using a higher order servovalve and uncompensated 
actuator models in the simulation are necessary. Possibly, the re-checking 
of the consistency of the numerical values used in the baseline simulation 
is warranted, although the discovery of a significant source of system 
damping seemed to make this unnecessary. Additionally, the investigation 
of linkage and load damping emphasized the need for continuing research 
into unsteady rudder system performance characteristics, including rudder 
lift and twisting moment variations with such variables as rudder turning 
rate, ship velocity, and ship turning rate. 
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V. CONCLUSIONS 



1. Analysis of unsteady rudder behavior identified a significant 
source of rate-dependent torque which had been neglected in the development 
of a submarine steering control system’s analog simulation. Programming 

of representative levels of feedback in the simulator virtually eliminated 
oscillations in rudder rate and load pressure drop which were thought to 
be unrealistic. 

2. The actual performance characteristics of a submarine steering 
control system under varying operational » conditions appeared to be loosely 
defined, and quantitative comparisons of system performance seemed to be 
necessary in order to continue the development of a realistic simulator. 

3. Studies showed that the inclusion of more complex models of the 
servovalve and actuator could have important effects on the flexibility 
and realism of the simulator, and that further work should be done to 
accomplish this end. 

4. Continuing theoretical and experimental research in the field of 
rudder behavior under unsteady conditions was suggested in order to improve 
the current concepts employed in steering system design and simulation. 
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Appendix A 

Explanation of Notation 



1 . Sat ( x 



■juL 

LL 



^ ->ul 
Sat(x) n 



ul 


x>ul 


X 


ll<x«ul 


11 


X<11 



This can be represented graphically as follows: 




Sgn(x) 



+ 1 x>0 

0 x=0 

-1 x<0 



This can be represented as follows: 
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